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ABSTRACT 



Beginning with linearized forms of the vorticity 
equation and the first law of thermodynamics applied to a 
simple baroclinic model including friction, an analytic 
solution is obtained which describes the time variation 
of the relative phase angle between the temperature and 
pressure waves in the atmosphere. It is found that for 
unstable waves the atmosphere tends toward a state where 
the temperature wave lags the pressure wave. The time 
variation of the amplitude and relative phase depends on 
the initial value of the phase difference, but the ultimate 
angle by which the temperature field lags the pressure 
field depends only on the drag coefficient, certain atmos- 
pheric parameters, and the wavelength of the waves. Figures 
are included to show the variation with time and initial 
phase difference of the amplitude and the relative phase. 
Also shown are the effects on the amplitude and terminal 
value of the phase difference due to variations in the drag 
coefficient and other parameters. 

The writer wishes to express his appreciation to 
Professor George J. Haltiner of the U. S. Naval Postgraduate 
School for his guidance, contributions and encouragement 
in this work. 
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1. Introduction 

In I960 Wiin-Nielsen [5J published the results of 
an investigation of the mechanisms in the atmosphere 
whereby the conversion of potential to kinetic energy 
is possible. He stated that for this conversion to occur 
it was necessary that the temperature wave lag behind the 
pressure wave. Using a simple baroclinic model with no 
friction he obtained an analytic solution for the relative 
phase between the temperature and pressure fields. Based 
on this solution he indicated how the relative phase angle 
changed with time for various initial values and how the 
early tendency of the pressure wave to amplify was affected 
by the initial value of the relative phase angle . 

The purpose of this work is to include friction in the 
simple baroclinic model and compare the results to the 
non-friction case. Also the values of the initial relative 
phase angle, wavelength, drag coefficient, Coriolis para- 
meter, thermal wind, and static stability were varied to 
investigate their effects on the amplification of the 
pressure wave, and the time variations including the limit- 
ing value of the relative phase between the temperature 
and pressure waves. 
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2. Development. 

The atmosphere is represented by a simple baroclinic 
model. The 750- and 250-mb levels are used as data levels, 
and the 500-mb level is taken as the level of interest. 

This gives a 4-layer model where CJ=0 at the upper boundary, 
P=0, (level zero); and Q=(^ , the frictionally induced value 
at level four, the lower boundary of the atmosphere. It is 
also assumed that the atmospheric motions are adiabatic. A 
zonal basic wind field which is independent of latitude and 
a linear function of pressure is utilized. 

The equations used are the vorticity equation in the 

form 

( 1 ) 



and an expanded form of the first law of thermodynamics for 
adiabatic motion given by Haltiner C2J 

$ + I2I 

Here c5“ is a measure of the static stability, 



e zp ip > 

and has been replaced by —\x/^P* It will be assumed 
that the basic flow pattern remains constant and that the 
perturbations grow, but no attempt will be made to account 
for the form of the necessary energy input for this to occur 
notwithstanding the dissipating effect of the friction. The 
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vorticity equation will be applied at 25C mb, level one, 
and at 750 mb, level three. The thermal equation is applied 
at 500 mb, level two. The wind at level i is denoted 
except that at level two the wind is taken to be 
The thermal wind for the layer of thickness AP centered 
about level two is given by 

For simplicity in the equations it was assumed that 
the frictional stress is directly proportional to the 
surface wind. The geostrophic wind approximation leads to 
the result Here C D is a drag coefficient, 

^ is geostrophic relative vorticity, g is gravity, and f 
is the Coriolis parameter. 

Applying the vorticity equation at levels one and 
three yields: 

and expressing ]y, and Q4 in terms of W and yields 
for level one the result 

+ (V -t-Vr) • VC^tlr ff) = A- YfiP J ( 3 ) 

and for level three 

+ ( 4 ) 

Here, the thermal equation is applied at level two to obtain 

~ (¥t +Vi‘Vk) / (e~ AP) 
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and the linear extrapolation i P use< i to obtain 

6)^. Adding equations (3) and (4) gives 

j?J +V-V(l+f) + W-V?T=-&&(T-2$. ( 5 ) 

Subtracting equation (4) from equation (3) produces 



ji% +Vr'V(^)-*V‘VJr 



+ ^(S-2?r). ( 6 ) 

To linearize these last two equations assume that K and 
\/ are given by 

V-(U+u)t+v7 W Vr-(U T *t^)t + 'lfS J 

where the U and U t represent the basic zonal flow which is 
a function of the pressure only, and the small u’s and v's 
represent the perturbation quantities which are independent 
of latitude, the y-axis. Also take h to be a basic part 
plus a perturbation, i.e., W-H(V) Now by neglecting 

terms involving products of perturbation quantities and 
replacing v by v=|^, where ^ is a stream function given 
by and “Vf- by Vr-'Q ^ t where is a stream function 

given by we obtain two third order linear partial 

differential equations in ^ and : 

jf " ^ i-iu, $ - 



Uere /?= ^/ K=-MC 0 jf9-fdP ) *«J S=ff~/23<r(6pf: 



( 7 ) 



Assume solutions of equations (7) and (8) of the form 

Y = yt e anc j 74- =-]4: a 

where 'X/y and r \j^ 0 are complex amplitudes >y/ £/=:^is the wave 
number for wavelength equal to L , and C is the phase 
speed. Substituting these two solutions into the differ- 
ential equations yields two homegeneous linear equations 
in ^ and ^ 




In order to have non-trivial solutions to this system of 
equations it is required that the determinant of the 
coefficients of and ys. 0 be equal to zero. Setting the 
determinant equal to zero and defining the 

Rossby wave speed, yields a quadratic equation in (C~Ci\) 
Solving this equation for (c-Cr ) gives 

C-Cr-F+lF -j# \jHi-i6 1 (11) 

G=LtUr(j£-ji?-3) 

A-o+$) 

By comparing equation (11) to the result given by 
Wiin-Nielsen /j5U it can be seen that the effect of including 
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friction in the model has been the addition of those terms 
involving A as a factor to the frequency equation. If 
the drag coefficient, C D , is set equal to zero, vanishes 
in equation (11) and the result is identical to Wiin- 
Nielsen's equation. 

Equation (11) gives C as a complex number where the real 
part determines the speed of propogation of the perturba- 
tions, and the imaginary part when multiplied by the wave 
number gives the amplification for the wave. Taking first 
the plus sign and then the minus sign with the radical gives 
two solutions which are designated as C+ and C„ respectively. 
The complete solutions for ''/'and are then given by: 

q £“( x-c.t) (12) 

y r = + Yro. (13 ) 

Next take as initial conditions at time t=0 

V'=//'e4'* 

where pf and A^, are real numbers and c«z, is a relative 
phase angle which is positive when the temperature field 
lags the pressure field. Then by considering equations (12) 
and (13) at time zero it follows immediately that 

By utilizing equation (9) for C+ and C_ we obtain two more 
equations for S ivin S a system of four 

equations and four unknowns. Solving the system gives 
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( 14 ) 



Vi^(KAre^-A'x-)/(x, -y.) Vfo + = -fy 0 + 
'V' 0 -=(£Ae ik -A / xl/(x.-Xj (15) 

where K-Ur-X ^£- j 

X*t ~ C-r~C R 4 A a.n.d X- — C- ~Cr "i- . 

Then by substituting the expressions for 'Y^ t > 

and yin- into equations (12) and (13) the complete solutions 
for Y'and ^ are obtained. 

In order to study unstable baroclinic waves in some 
detail, it is desirable to have the expressions for ^ and 
'^ 7 * i n slightly different forms. Since C + and C_ are complex 
they can be expressed as C + =Y+iW and C_=M-fiN. Then the 
expression for ^ becomes 

y / - yX + e ^ e ^cx-xtA) 

Now define B+ , B- , <X+ , and by 

Tfce**-ae'* ' 1S1 

Then 

y = + /2~e L<k -e. 

which after some manipulation can be written as 

y = /3 q iC/ix-MG \+F)t t&j 

Here j&is given by 

B’ 1 '= Ql +B- ■+ 2 [(&+-<£.) -■?]> (17) 
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(IS) 



<£> is given by 

_ /3+ sA ( -£J -+ 13- 0^-^) . 

G * Q-ctt (&+%) J 
is defined as ir=£fp6«$t where f = sf^f+Jp and 
£=arctan (G/H); F, G, H, and A are defined immediately 
following equation (11). The treatment for 0^-is similar 
to the above procedure for ")4i Define Qr+ l Cr- y ^ 

-J K.e'n&.e*'* 

and then obtain 

/4>x-//rc„ vv» * + <£tJ 

with ^.given by 

Hr ^ f 2 f3n Hr- &*■!&•+ "* 5r-)"" *2^7, ( 2 °) 

and c^y given by 

/- ~ -Hr* + !3 t- &*.(*&- + %) . (2i) 

We now have obtained expressions which give the variations 
of O^and ^as simple sinusoidal waves having amplitudes and 
phase angles which are functions of time. 

We now want to determine how the relative phase angle 
between the thermal and pressure waves, { ) changes 

with time. Considering the exponential nature of J3^ f3r+> 

and /3y_, it can be shown from equations (IS) and (21) that 
as time approaches infinity, ^approaches (V^--^") and <p> r 
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approaches (<J^. + — 3r). Therefore, in the limiting case 

— <$>) =(5^. — ^). Now by suitable combining equations 
(14), (15), (16), and (19), it can be shown that 

_ Ut O* -r +7\P^-fr & 1 ) +2 A 

Equation (22) was obtained by considering only the 
amplifying part of the wave. Unstable baroclinic waves 
tend toward a state such that the limiting phase angle 
between the thermal and pressure waves, (^-y. — 5^. ) , is given 

by $ =FMCT/Ofi/(Vr,l, e-jt, Co, 3) 

where the function is given by the right side of equation 
(22). It is to be noted that the initial value of the 
relative phase difference,^, does not appear in the 
expression for the limiting value of ( cpr — 9^ ) . 
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3. Discussion of Results. 

Equation (22) was solved numerically and the results 
are plotted in figure 1, which gives ( <£?+ — <p + ) as a 
function of wavelength and drag coefficient » Using a 
surface geostrophic wind of 10 m/sec, the values of C D in 
figure 1 correspond to frictional stresses of from zero to 
20 dynes/cm^. 0^=0.004 dyne sec/cm^ corresponds to the 
value given by Sutton f43 for grass about 15 centimeters in 
height. The value of the static stability parameter;, , 

of 2.90x10“^ cmVdyne 2 corresponds to the value given by 
Gates f 13 as the mean value over the United States in 
January. Those curves in figure 1 for which *$).) is 

zero when the drag coefficient is zero, are the curves for 
which term H in equation (11) is positive, i.e. these 
wavelengths are the stable waves since the requirement for 
instability in the frictionless case is that H be negative. 
The apparent amplification depicted in figure 2 for what 
should be stable waves with C & =0 can be explained by 
reference to equations (16) and (17) from which it can 
be seen that when C + and C_ are real numbers, the maximum 
possible value of B is given by b 2=B+ + B 2 + 2B + B_. The 
time variation of B for two stable waves, 3,000 km and 
12,000 km, with C D =0 is shown in figure 3. As is apparent 
from equation (17), B is a sinusoidal function of time, 
oscillating about the initial value. Figure 2 shows the 
effects of variations in the drag coefficient on the value 
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of the amplitude for the wavelength band from 2,000 km to 
14,000 km. The values of B used to construct the figures 
are the values attained at the end of the first 36 hours 
for those wavelengths having {<P T + — <P+ )>0, for 0^=0, in 
figure 1. For the stable waves, due to the oscillatory- 
nature of B mentioned above , the maximum value of B 
attained during the 36-hour period was used. For each wave- 
length the value of the initial phase difference between 
the temperature and pressure fields,^, was such as to 
maximize the amplification in 36 hours (oc given to the 
nearest 45 degrees). As is obvious from figure 2, the 
effect of increasing the drag coefficient is to reduce the 
amplification for all wavelengths considered. This effect 
of friction is in agreement with the results presented by 
Holopainen In comparing the results it should be 

noted that the ordinate in figure 2 is wavelength, whereas 
the ordinate in Holopainen' s figure is a function of wave- 
length, static stability, and latitude. In figure 2 it 
can be seen that the wavelength for which the amplifi- 
cation is maximum decreases as the drag coefficient 
increases. 

Figures 4, 5, and 6 depict the effects of variations 
in thermal wind, static stability, and latitude on the 
maximum value of the amplitude for the same wavelengths as 
were used in figures 1 and 2, with C D =0.004 dyne sec/ cm?. 

The central curves in figures 4, 5, and 6 are identical and 
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are the same as the middle curve in figure 2„ From figure 4 
it is apparent that as the thermal wind increases the 
amplitude increases for all wavelengths, and that the wave- 
length for which the amplitude is a maximum also increases . 
Figure 5 shows that the effect of increasing the static 
stability is to decrease the amount of amplification for all 
wavelengths, except that for wavelengths greater than 
about 13,000 km the amplitude increases slightly. Also as <5~ 
increases the wavelength of maximum amplification increases. 
The effects of moving the latitude northward are, as 
depicted in figure 6, to increase the amplification of the 
entire waveband, and to shift the wavelength of maximum 
amplitude toward lower values. It appears from the figures 
that increasing the static stability and decreasing the 
latitude have essentially the same effects on the amplifi- 
cation, and that decreasing the thermal wind corresponds 
to increasing the drag coefficient in the resulting effect 
on the wave amplitude. One might conclude from the 
information presented in the figures that at low latitudes 
with a high value of the static stability parameter and a 
low thermal wind the pressure waves would not tend to 
develop, and that at high latitudes with low static 
stability and high thermal wind the tendency would be toward 
amplification of the pressure waves. 

Wiin-Nielsen discusses the validity of applying the 
results of a linearized treatment of the dynamical equations 
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to obtain a limiting case for very long times. He justified 
the application on the basis that the initial trends were 
toward the situation given as the limiting case. With 
similar logic in mind, equations (17), (IS), (20), and (21) 
were used to numerically investigate the values of the 
amplitude factors B and B r , and the relative phase difference 
(<p r -cp) for periods of time of less than 72 hours. It 
should be noted that while the limiting value of («j^--<^>) is 
independent of the initial phase difference , oc, the short 
period variations of B, B T , and are very much 

dependent on the value of oc. Figures 7 and £ were 
constructed for a 5,000 km wavelength wave to depict the time 
variation of (<^--<£) and B, to show the dependence on of 
the time variation, and to show the fact that the terminal 
value of (^ — ^) is independent of the value of <x. „ 

Because of the very similar behavior in the variations 
of B and B r with time, only B is shown in the figures. 

Figures 7 and S may be compared with similar figures 
presented by Wiin-Nielsen £53 , using the same wavelength, 
for the frictionless case. In comparing the two sets of 
curves consideration should be given to the differences 
in the values of the atmospheric parameters used in the 
calculations performed to construct the curves in figures 
7 and S, and the values used by Wiin-Nielsen. One 
conclusion that can be drawn from the figures is that 
increasing the drag coefficient decreases the time required 
for (cp T -<P ) to approach its limiting value. 
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The number at the top of each curve is the wavelength xlO 
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Figure 




Figure 2 
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La.t/ tnJe = 4*5/V Ur - 20 m/se. c. 2,70*10 
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(^Pt — C P) decrees 




T~/ <2 Hours. 

Figure 7 
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